function [CE_HHW_man, share_floor] = CE_HHW_crra(Wi,iX,MU,CAPH,Nstates,rho,T,tolFP,maxreps,sigma,floor)
    Y1 = Wi(:,1);
    Cg1 = HHW_cg(Y1,MU,CAPH,Nstates,rho,T,tolFP,maxreps);
    C = Consumption_HHW(Cg1,iX);
    share_floor= mean(mean(C<floor));
    C(C<floor)=floor;
    u = 1/(1-sigma)*(1000*C).^(1-sigma);
    Num = mean(sum(u.*(iX~=(Nstates+1)).*power(rho,(1:T)-1)',1));
    Den = mean(sum((iX~=(Nstates+1)).*power(rho,(1:T)-1)',1));
    CE_HHW_man = ((1-sigma)*(Num/Den))^(1/(1-sigma));
end